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ABSTRACT 

We present a systematic evaluation of the agreement between the observed 
radii of 90 well-characterized transiting extrasolar giant planets and their cor- 
responding model radii. Our model radii are drawn from previously published 
calculations of core-less giant planets that have attained their asymptotic radii, 
and which have been tabulated for a range of planet masses and equilibrium 
temperatures. (We report a two-dimensional polynomial fitting function that 
accurately represents the models). As expected, the model radii provide a sta- 
tistically significant improvement over a null hypothesis that the sizes of giant 
planets are completely independent of mass and effective temperature. As is well 
known, however, fiducial models provide an insufficient explanation; the plane- 
tary radius anomalies, TZ = -Robs ~ -Rpred, are strongly correlated with planetary 
equilibrium temperature. We find that the radius anomalies have a best-fit de- 
pendence, TZ oc Tgg, with a = 1.4 ± 0.6. Incorporating this relation into the 
model radii leads to substantially less scatter in the radius correlation. The 
extra temperature dependence represents an important constraint on theoreti- 
cal models for Hot Jupiters. Using simple scaling arguments, we find support 
for the hypothesis of Batygin and Stevenson (2010) that this correlation can be 
attributed to a planetary heating mechanism that is mediated by magnetohy- 
drodynamic coupling between the planetary magnetic field and near-surface flow 
that is accompanied by ohmic dissipation at adiabatic depth. Additionally, we 
find that the temperature dependence is likely too strong to admit kinetic heat- 
ing as the primary source of anomalous energy generation within the majority of 
the observed transiting planets. 
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Introduction 



The disparate radii of the transiting extrasolar planets have defied a straightforward 
explanation. With nearly 100 well-characterized examples now known, it is clear that where 
Jovian planets are concerned, variations in mass and stellar insolation are responsible only 
for a fraction of the observed range in planetary sizes. 

Radius a nomalies emerged Avith th e discovery of the first transiting extrasolar planet, 
HP 209458b (ICharbonneau et al.1 120001). Structura l evolu tionary models computed by, e.g., 
Bodenheimer et al.l (120011 ) and iGuillot fc ShowmanI (120021 ) suggested that a H-He dominated 
planet with HD 209458b's mass, insolation, and age should have a radius of order R ~ 
1. 1-R.inp. This figu re is startlingly at odds with the observed value, R = 1.38 ± 0.02i?jup 



(ISouthworthll2010f ). 



A decade of additional discoveries has made it clear that HD 209458b 's radius anomaly is 
by no meai is anomalous, and the st atus of discussion ir i the field i s well- covered i n the recent 
reviews by IShowman et al.l (120101 ) , iBurrows fc Ortoru (120091 ) , and iBaraffe et al.l (120101 ) . The 
extraordinary variation in observed radii in the planetary population is indicated by Figure 
1. This diagram charts the measured radii and uncertainties for the 90 transiting planets 
with accurately measured masses against the planets' orbit-averaged effective temperatures, 
given by 

2a / fl 



Teff = ( - 



:ii 



where i?^ is the stellar radius, a is the semi-major axis, e is the orbital eccentricity, is 
the stellar effective temperature, and zero albedos have been assumed. Only planets with 
O.lMjup < Mpi < 10.0 Mjup were used£| 

It is evident that the observationally well-characterized transiting planets exhibit a wide 
range of radii for a given mass, and a number of models have been ad vanced to e xplain this 
dispersion. Planets that are smaller than expected, with HD 149026b ( 1Satdl2005l ) providing 
the canonical example, have generally had their small radii attributed to a large fraction 
of heavy elements in their interiors. For planets such as HD 209458b, which are larger 
(and sometimes much larger) than expected, explanatory models must generally appeal to a 
cryptic source of heating at adiabatic depth. A number of heating mechan isms have been pro- 



posed , including dissipation stemming from tidal orbital circularization (iBodenheimer et al. 
200 ll), ' 'kinetic heating" in wh ich wind energy is converte d into heat (IGuillot fc ShowmanI 



2002i) . and Ohmic dissipation (IBatygin fc Stevensonll2010l ). This letter shows that current 



^AU observational data taken from www.exoplanet.eu, accessed 12/03/2010. 
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observational evidence indicates that the radius anomaly, the difference between the observed 
radius and that predicted by theoretical models, displays a well-defined correlation with tem- 
perature (from equation [1]). Further, this temperature dependence can be interpreted to 
favor a magnetohydrodynamic mechanism in which Ohmic heating and attendant magnetic 
braking of the velocity field at sufficient atmospheric ionization fraction play significant roles. 
Our approach relies on scaling arguments, as we believe that the extant observations are not 
yet sufficient to permit the discrimination between detailed models containing multiple free 
parameters. 



2. A Comparison Between Structural Models and the Observations 

Naively, one might expect that the radius of a mature gas-giant planet is primarily deter- 
mined by its mass and by the amount of radiative energy that it receives from its star. This 
conjecture can be tested by evaluating model radii, Rm^, of "baseline" evolutionary models 
for H-He composition planets spanning a range of masses and insolation, and comparing with 
those of known corresponding transiting planets, i?o,- If the model has explanatory power for 
the aggregate of known planets, then it should produce a statistically significant decrease 
in the quantity. 



2 _ 1 \^ {Rrrii Roi^^ 

m 



in comparison to Xnuih obtained by replacing 1?^^ with Rav = 1.2i?jup, the average observed 
radius for transiting planets having O.lMjup < Mp\ < lO.OMjup. In the above equation, 
Nf = 2 is the number of free parameters (Mpi, Teg) in the explanatory model. 



Our baseline models were published by iBodenheimer et al.l (120031 ): hereafter BLL. As 



described in BLL, radius estimates, Rrm, were computed with a Henyey-type planetary 
structure calculation, and the reader is referred to that paper for details regarding the input 
physics and assumptions. The model radii were tabulated at 4.5 Gyr for a grid of assumed 
Mpi and T^s- We ignore the small dependence of radius on age for mature planets and use 
bilinear interpolation to obtain an estimate for R^i at given Mpi and T^s. Our estimates are 
drawn from BLL's model sequence of core-free solar-composition planets with no anomalous 
energy sources. 

For the 90 transiting planets, we find Xm = 23.5 and Xnuu — 32.6. Not surprisingly, this 
result indicates that the baseline structural models can explain some, but by no means all, 
of the the observed variation in planetary radii. As an alternative to bilinear interpolation 
between table values, it can be useful to have a simple fitting relation, Rpi{Mp\, Tes)- Defining 
m = logiQ(Mpi/Mjup), and t = Tgff/lOOO, we find that the two-dimensional polynomial fitting 
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function, 



Rp\/Rinp = 1.08417 + 0.0940857m - 0.242831 + 0.0947349 
+0.0387851 1 + 0.00243981 mt - 0.0244656 niH + 0.0130659 mH 
+0.0240409 - 0.0419296 mt^ + 0.00693348 m V + 0.00302157 m V , 



(3) 



provides an acceptable approximation to the BLL baseline structural models throughout the 
region where O.lMj^p < Mpi < lO.OMj^p, and 100 < T^s < 25000 

For each planet, we define a radius anomaly, TZi = Rm, — Roi, and look for correlations 
between 71 and other measured quantities (such as Tes, T*, M^, etc). Many authors, (e.g., 
Enoch et al. 2010) have noticed that planetary radii tend to swell dramatically with in- 
creasing insolation. Figure 2 illustrates the s ignificant correlat ion between TZ and planetary 
Teff. Using a bootstrap replacement method (jPress et al.lll992[ l. we find a best fit power-law 
dependence, 

n oc T^^^'-' . (4) 



Among the various mechanism that have been invoked to explai n the radius anoma- 
l ies, we expect that both Ohmic heating (iBatygin fc Stevenson! |2010[ ) and kinetic heating 
( IGuillot fc Showman! |2002[ ) should show a positive correlation between TZ and Tes- We can 
ask, furthermore, whether the measured exponent, a = 1.4 + 0.6, is consistent with either of 
these proposed mechanisms. 



In the treatment of !Batygin fc StevensonI ( !2010! ). energy deposition in the planet is 
approximated by integrating over the resistivity in each mass element 



E 




-dV . 



(5) 



where J is the current density and cTr(r) is the conductivity in the near-surface (but still 
adiabatic) layers where the resistive heating occurs. For order-of-magnitude purposes, they 
ignore the radial dependence, ar{r) ~ 0"^,. Further, it is assumed that the relevant currents 
are induced by the planet's intrinsic dipole field, B, so that 



J = cr,.(v X B) 



and 



E oc (xJvPlBp, 



(6) 



We want to estimate the temperature dependence, £'(Tefj), implied by the Batygin-Stevenson 
model. First, we need an estimate of how wind kinetic energy at adiabatic depth scales with 



^Applying this polynomial relation in place of bilinear interpolation yields Xm' — 18.30 for the aggregate 
of 90 transiting planets. 
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temperature. One possibility is that there is a constant partitioning between thermal en- 
er gy and kinetic energy, an a ssumption that conforms with the kinetic heating hypothesis 
of iGuillot &: ShowmanI (120021 ) . which posits that a fraction, rj, of the total flux received by 
a planet is converted into kinetic energy by atmospheric pressure gradients, and that this 
energy is dissipated at depth by Kelvin-Helmholtz instabilities. In this case, one has v ocT^ 



eff- 



Alternately, one can examine three-dimensional hydrodynamical models for irradiated giant 
planets that experience a range of insolations. These models suggest a much weaker depen- 
dence of wind speed on T^s at large optical depths where energy diss ipation can effectively 
couple to the planetary structure. For example, the lLewis et al.l (120101 ) model for Gliese 436b 



(Mpi = 24M0, Tcfj = 707 K) which adopts a solar metallicity atmosphere, finds v ~ 100 ms~^ 
at P = 10 bar depth. This wind velocity is ne arly identical to the average wind velocity at 
10 bar depth found by IShowman et al.l (120091 ) for a solar metallicity atmospheric model of 



HD 189733 (M = 1.2Mjup, = 1202 iT). For our rough analysis, we choose v oc T^^'^, 
bearing in mind that the temperature's power-law index might reasonably fall anywhere in 
the range from to 2. 



Second, we follow iBatygin fc StevensonI (120101 ) and assume that the Elsasser Number, 
A = (TjB^/2pf2, is roughly constant across the aggregate of transiting planets. (The quantity 
(Tj is the conductivity appropriate to the interior regions where the B field is assumed to 
be generated by dynamo action, p is the density of the planet, and Q is its angular spin 
frequency). To the precision of our analysis, this relation holds for the bounding cases of 
Jupiter and the Sun; Jupiter has a magnetic field roughly ten times the solar value, and 
has a spin period that is roughly 100 times shorter. Both bodies have similar densities and 
are highly conductive in their interiors. With A and (jj assumed constant, we can remove 
the magnetic field dependence, yielding E oc arpfiTes ■ If we further assume synchronous 
rotation (which should hold for hot Jupiters on circular orbits), then l/f2 cx a^/^ and a ocT^ 
(via equation [1]), giving E oc cr^Tjj. 



eff 



We need to identify how the dissipation-layer conductivity, 0"^, depends on temperature. 
The conductivity is directly tied to the ionization fraction, which in turn is determined by 
the Saha equation 

27rh^ J 



n- He 



9j 



■exp [—Ij/kT] 



(7) 



where the j chemical species typically include Na, K, Li, Rb, Fe, Cs, and Ca, and T is 
the local temperature. The quantity Ij is the ionization potential (for the first electron to 
be ionized) and ranges from 4.3 eV (for K) to 7.9 eV (for Fe), whereas g^, and gj are the 
degeneracies of states for the ions and neutrals, respectively, of species j. Note that ng is 
the total electron number density (since all electrons can interact with each species of ion). 
In the relevant temperature regime, ionization fractions are low, so that n'j ^ rij. Summing 
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over all species and simplifying yields 



Up 



m, 



,kT\ 



3/4 



27rn2 



1/2 



(8) 



where fj = 2g^/gj. 



Next, we note that the electrical conductivity, 0"^, appropriate to the dissipative region 
has the form 



Or 



n rripA \8kT 



\8kT) 



(9) 



where A is a weighted cross section, see, e.g. iTipler fc Llewellynl (120021 ). The quantity of 
interest is the index p defined by 



P 



T d(Tr 



1 T dup 

2 ^ 



(10) 



ar dT 

The index p allows us to estimate the temperature dependence of the conductivity, ar- Using 
our expression for ng, we obtain 

-I -1 



1 T 
^=4+2 



^/jnjexp[-/,/A;r] 



d_ 
df 



^fjUjexp [-Ij/kT] 
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After differentiation, this equation can be written in the (apparently) simpler form 



1 



P = 



{I) 



4 2kT 



:i2i 



where we have defined 



{I) 



^/^n.exp [-Ij/kT] 



^IjfjUjexpl-I^/kT] 



(13) 



Since kT -C Ij, the exponential suppression factors dominate the behavior of this function. 
As a result, we can consider the largest term to dominate (which will be the term with the 
smallest ionization potential /*): 

1 T 

(14) 



1 

p ~ - 



4 2kT 

For T = Teff = 2000 K, the quantity 2kT = 0.345 eV. If = 4 eV, then p ^ 12. It's 
clear that i?(Tefj) depends sensitively on Tgfj, and that for a given choice of the index p 
increases sharply with T^s (for /* = 4 eV, p=15.7 at 1500 K, and p=9.5 at 2500 K). 

A similar result is obtained by Batygin & Stevenson (2010), who present a simple but 
accurate model for the conductivity as a function of radius, where a oc exp[— (r — R)/H], 
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where H is the conductivity scale height (which is about twice the thermal scale height). 
Since H ^ T, the index p ^ {r — R)/H, which is of order 10. 

Next, we need to establish the dependence of the radius anomaly, TZ on the internally 
generated power, E. This depends on the planetary mass and the size of the assumed 
core - it's easier to inflate a low-mass planet without a core. For a typical hot Jupiter, 
such as HAT-P-13, models based on the BLL prescription indicate that by increasing E 



more than tenfold, from 3.5 x 10^ ^ ergs 



^ to 5.6 X 10^^ ergs ^, one inflates the radius from 



R = 1.2i?jup to i? = 1.36i?jup (IBatygin et al.ll2009l ). T his single examp le suggests that 
71 depends on E raised to a small fractional power, x. iGu et al.l (12004 ) have computed 
sequences of evolutionary models (again using a method similar to that described in BLL) 
which systematically explore how planetary radii depend on internally generated power. We 
obtained a fit to their sequences and find that an index a; ~ 1/6 applies over a wide range of 
masses for planets having log(i?pi/i?jup) < 0.3. This regime covers nearly all of the known 
transiting cases. Using p=l2, we arrive at a predicted dependence 

, (12+4)/6 



7^ oc T" ^ T 



eflf 



m 2.666 



(15) 



for which the index a is more that one power of temperature steeper than suggested by the 
currently observed aggregate of transiting planets (JZ oc T^-^). 

Caution regarding the quantitative aspect of the dependence given by equation (15) is 
in order. Order-of-magnitude scaling arguments can at best give only rough insight into the 
physical balance that determines the radii of the hot Jupiters. As mentioned above, there 
are considerable uncertainties associated with the wind speeds at adiabatic depth, and the 
treatment leading to equation (14) is highly idealized. We have assumed, furthermore, that 
the Ohmic heating rate is directly proportional to the conductivity, ar, in the weather layer. 
If the ionization fraction is large enough, a more fully self-consistent treatment must account 
for the back-reaction onto the wind velocity, v, generated by the Lorentz force 



civ 1 , -r-. 
p— = -(Tr(v X B) X B . 
at c 



(16) 



Indeed, iPerna et al.l (l2010bl ) and iPerna et al.l (l2010al ) have suggested that the bulk 
Lorentz force acts as a drag term on the velocity field. In this event, higher tempera- 
tures increase ionization, generating Lorentz forces that brake the velocity field. Reduced 
wind-speeds lower the rate of Ohmic dissipation, which in turn would decrease the tem- 
perature power-law index a. A complete understanding of this process requires a full 3D 
MHD approach. While such simulations are not beyond current computational capabilities, 
the lack of detailed knowledge of the atmospheric conditions likely renders them somewhat 
premature at this time. 
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3. Discussion 

Transiting extrasolar planets tend to liave radii tliat are larger tlian tliose predicted 
by structural models parameterized by Mpi and Teg. The resulting radius anomalies, TZ, 
are strongly correlated with effective temperature, with TZ oc T^^^ providing a maximum 
reduction in variance. The simple scaling arguments of the previous section suggest that 
Ohmic heating may play a significant role, and that further progress may be made as self- 
consistent MHD-mediated weather models are constructed. 

A very rapid increase in E with temperature pos es difficulties for some pos ited ex- 



planations for the radius anomalies. For example, the iGuillot fc ShowmanI (120021 ) kinetic 
heating mechanism proposes that a constant fraction, rj ~ 0.01, of the total fiux received by 
a planet is converted into kinetic energy by atmospheric gradients, and that this energy is 
dissipated at adiabatic depth. While the flow patterns on the surfaces of strongly irradiated 
extrasolar planets are still a matter of debate, a E oc Tjj scaling is implied by the kinetic 
heating hypothesis, leading to 7^ oc T^/^, which (to > la confldence) is weaker than the 
observed dependence. This argument does not imply that kinetic heating is absent, rather, 
the implication is that it is unlikely to be the primary mechanism that determines the radius 
anomalies. 

Finally, we note that substantial variance in the observed radii remains after the TZ oc 
Tgg"^ scaling has been removed. Inclusion of the best-flt (3T°^ dependence into the model 
radii increases Nj from 2 to 4, while causing Xm measured with Equation [2] to decrease 
from Xm = 23.5 to Xm' = 14.8. Clearly, processes other than Ohmic heating are also 
contributing. For example, our base-line models unrealistically assume solar-composition 
planets. Both Jupiter and Saturn are substantially super-solar, and the small radii of planets 
such as HD 149026b and CoRoT 8b can only be understood if elements heavier than H/He 
contribute more than 50% of the planetary mass. Figure 3 indicates that there is evidence for 
a signiflcant correlation between the residual radius anomalies and host star metallicities. 
This correlation is generally interpreted as evidence that metal-rich protoplanetary disks 
lead to planets with larger core masses. Furthermore, a fraction of the known transiting 
planets have eccentric orbits, meaning that interior heating from tidal orbital circularization 
must be ongoing in some of the observed planets. For modest eccentric ities, and at a given 



time- average tidal quality factor, Q, tidal heating rates scale with ( iMurray &: Dermott 



19991 ). This strong dependence could be a determining factor in producing the residual 
radius anomalies for many of the observed planets. Furthermore, even if a planet's orbital 
eccentricity is currently zero, it may still be inflate d as a consequence of th ermal inertia from 



a tidal heating episode that occurred in the past ( llbgui fc Burrows! 120091 ). 



The number of transiting planets with well-determined properties will increase rapidly 
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in the coming years. An expanded catalog of planets, in conjunction with improvements 
to the structural models, will allow empirical quantities, such as the power-law index "a" 
that we've used here, to be determined with increased confidence. In addition, follow-up 
campaigns from platforms such as JWST will improve our understanding of the atmospheric 
conditions on short-period planets, and we can look forward to a time when the construction 
of detailed, fully self-consistent 3D MHD climate models are a sensible response to the quality 
and depth of the observational data sets. 

The authors acknowledge useful discussions with Konstantin Batygin, and we are grate- 
ful to an anonymous referee for a prompt and insightful report. This work was funded by 
NASA grant NNX08AY38A and NASA/Spitzer/JPL grant 1368434 to GL, and by NASA 
grant NNX07AP17G to FCA. 
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Fig. 1. — Radii and estimated effective temperatures for 90 transiting extrasolar planets 
with well-determined masses and radii. Circle size is in proportion to planetary mass, and 
the color of the inner circle indicates metallicity. 
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Fig. 2. — Radius anomaly, TZ, vs. planetary effective temperature, Tgfj, (see Eqn [1]) for 
90 transiting extrasolar planets. The red line charts the best error-weighted power-law fit 
to the data (JZ oc T", with a = 1.4). The 10,000 light gray lines show analogous best-fit 
power-laws to bootstrapped data sets in which the contributing planets are redrawn with 
replacement from the original data. The 5**^, 15**^, 85*^^ and 95**^ percentile bootstrap fits are 
indicated with black lines. We have adopted the 15*^^ and 85*^ percentile fits as an estimate 
of the 1 a confidence limits on a, where we find 6a ~ 0.6. For each planet plotted, the size 
of the associated variable circle is proportional to planetary mass. The gray-scale color of 
the inner circle indicates host-star metallicity. 
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Fig. 3. — Residual radius anomaly vs. host star [Fe/H] for 90 transiting extrasolar planets. 
The best-fit power-law and bootstrap realizations are as in Figure 2. 



